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Abstract. We report on 3.6 cm VLA and VLBA obser- 
vations of YZ CMi and AD Leo, two nearby dMe stars. We 
resolve YZ CMi and can fit a circular symmetrical gaus- 
sian component of FWHP of 0.98 ±0.2 mas, correspond- 
ing to an extent of the corona above the photosphere of 
1.77 x 10 10 ± 8.8 x 10 9 cm or 0.7 ±0.3 R+ (R* refers to 
the photospheric radius). We obtain an estimate of the 
brightness temperature of 7.3 x 10 7 K, which is consis- 
tent with that expected from gyrosynchrotron emission. 
For AD Leo the emitting region is unresolved. We there- 
fore set a conservative upper limit to its diameter of 1.8 
times the photosphere diameter, which leads to an extent 
of the corona above the photosphere of < 2.8 x 10 10 cm or 
< 0.8 We compare the radio emitting dMe stars with 
measured sizes with the Sun and conclude that these ac- 
tive stars have much more extended coronal radio emission 
than the Sun. The VLBA position of YZ CMi has been 
found to differ by 32 mas from the positions calculated 
from the Hipparcos catalogue. The discrepancy is caused 
by large errors in the listed proper motion. An improved 
value is given. 

Key words: Radio continuum: stars-stars: coronae- 
Astrometry-stars: individual:YZ CMi, AD Leo 



1. Introduction 

Hot coronae are the enigmatic link between cool, convec- 
tive stars and their environment. While most of this hot 
plasma (> 10 6 K) is contained by magnetic fields, some 
may escape due to overpressure as a hot stellar wind along 
open field lines. The extent of intense coronal emissions in 
X-rays or radio waves thus yields a lower limit to the size 
of closed magnetic loops, i.e. the size of the stellar mag- 
netosphere. 

Radio emission from solar active regions at 3 cm wave- 
length originates from altitudes above the photosphere 



between 5 to 10 x 10 8 cm, or 0.006 - 0.015 R Q , (e.g. 



Aschwanden et al. 1995). It is generally attributed to ther- 



mal gyroresonance emission. The size of the radio image 
of the Sun increases at longer wavelengths. In soft X-rays 
and EUV line emissions, coronal loops reach 3 x 10 10 cm 



(0.5 Rq) (e.g. Bturrock et al. 1996). Magnetic loops in ex- 
cess of 1.4 x 10 11 cm (2 Rq) have been observed in metric 
U bursts (Labrum & Stewart 1970). Coronograph obser- 
vations in white light, however, do not regularly see loops 
of that size, suggesting that they are shortlived. In gen- 
eral, large magnetic loops are confined to low latitudes 
and, consequently, the radio image of the Sun is more ex- 
tended in the equatorial direction. 

Very Long Baseline Interferometry (VLBI) has made 
it possible to study the size and shape of main-sequence 



stellar coronae. At 3.6 cm, Benz ct al. (1998) have resolved 
the dMe star UV Cet B at an higher and variable flux level 
into two components separated by 4.4 x 10 10 cm (4.4 R+). 
One of the components was spatially resolved with the 
size of about the stellar photosphere. The orientation of 
the two sources was found to lie along the probable axis 
of rotation, strongly suggestive of coronal enhancements 
extending at least 2.1 x 10 10 cm (2.1 R+) above the poles. 
Thus, UV Cet B differs considerably from the Sun in size 
and shape of the radio corona. 

More VLBI observations of dMe stars have been re- 



ported at 18 cm wavelength. For YZ CMi, Benz fc Alel 



(1991) found an upper limit of 8.7 x 10 u cm, suggesting 
an extent of a circular corona above the photosphere of 
less than 1.8 x 10 10 cm (0.74 R+). Two size measurements 
of AD Leo by Benz et al. (1995) suggest a coronal extent 
of less than 2 x 10 10 cm (1.0 R+) and 5.4 x 10 10 cm (2.6 
R+) above the photosphere, assuming a circular, concen- 
tric shape. A small total size of less than 4.9 x 10 10 cm for 



EQ Peg (Benz et al. 1995) refers to an observation during 
a totally polarized flare, and an extremely large size of 
2.2 x 10 11 cm was reported for the dMe close binary YY 
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Gem ( |Alef et al. 1997Q . 

The radio coronae of active, rapidly rotating dMe stars 
have been modeled to determine the emission process. It 
is generally agreed flGiidel 1994j ; |Whitc ct al. 1989Q that 
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weakly polarized radio emission is produced by the gy- 
rosynchrotron mechanism of a population of mildly rel- 
ativistic electrons. This mechanism, however, cannot ac- 
count for polarizations exceeding about 50% which may 
originate from a coherent process. Circular polarizations 
of up to 80% have been reported for YZ CMi at 20 cm 



and 18 cm ( [Lang fc Wilson 1986} |Bcnz fc Alcfl 1991). 

Here we report on the results of VLBA experiments of 
YZ CMi and AD Leo at 3.6 cm. These are well known, 
nearby, young radio stars close to zero-age main sequence 
(ZAMS). Some of their general properties are listed in 
Table [TJ. 

2. Observations and data reduction 

The single dMe stars YZ CMi and AD Leo were observed 
on 1997 April 18/19 and 1996 December 12, respectively, 
with the VLBA and the phased-up VLA as a joint VLBI 
system yielding an angular resolution of better than one 
milliarcsecond (mas) . 

The VLA was also available in its normal interferom- 
eter mode and was used for total flux measurements at 
very small baseline length which allowed us to monitor 
the changes in total flux density and polarisation of both 
stars through the observations. 3C286 was used for the 
flux calibration of both stars. 

The VLBI observations used a bandwidth of 8 MHz 
in both left and right circular polarization at 8.41 GHz 
(3.6 cm) and two bit sampling, giving a data rate of 
128 Mbit/s. In order to reliably image such weak sources 
( < 3.0 mJy) the phase referencing technique was used 
(Beasley & Conway 1995) in which we switched between 
the target and a bright calibrator (quasar) in cycles of 
two to three minutes. For our targets YZ CMi and AD 
Leo these calibrators were 0736+017 and 1022+194 re- 
spectively (see Table ||), at 2.1 and 1.4 degrees separation 
from the target. 

The amplitude calibration was performed with AIPS, 
followed by editing in DIFMAP (in order to flag more 
precisely single, bad visibilities of the VLBA data) and 
by a continued analysis in AIPS (including all the map- 
ping) . We first made hybrid maps of our calibrators using 
closure phase methods in order to determine their struc- 
ture. We then determined the atmospheric contributions 
to the phase of the calibrator data. With these solutions 
it was possible to phase-correct the stellar data. We then 
made preliminary wide maps in order to find the stars and 
phase-rotated the data in order to bring the target stars to 
the image phase center. Finally we deconvolved the stellar 
images using CLEAN. 

As a check of the reliability of the phase calibration, 
our observations also included a second bright calibrator 
for each star (0743-006 for YZ CMi, and 1013+208 for 
AD Leo, see Table |2|). The maps of these secondary cal- 
ibrators, made using the phase solutions found towards 
the primary calibrators showed a dynamic range of about 



20:1. Since the separation between the two calibrators is 
approximately the same as that between the star and each 
calibrator, we can conclude that our stellar images have 
also a 20:1 dynamic range and hence are in practice noise 
limited rather than dynamic range limited. We also pro- 
duced selfcalibration maps and datasets of the secondary 
calibrators. The plots of the fringe amplitude against u,v- 
distance show the same fall off than the ones produced 
from the phase referenced datasets. This confirms the 
value of the dynamic range of the images of the targets. 

The source centroid positions of the phase referenced 
maps of the secondary calibrators were within 0.5 and 
0.2 mas of the correlated positions (a-priori positions) for 
0743-006 and 1013+208 respectively (see Table |). This is 
consistent with the claimed accuracy of those correlated 



positions (0.5 mas, Johnston et al. 1995). This test gives 
us confidence in the astrometric accuracy of our stellar 
positions (see Sect. 3.4). 

3. Results 

3.1. Total flux monitoring 

YZ CMi was found at a surprisingly high, slowly decreas- 
ing flux density (average 2.9 mJy) during the whole ob- 
servation (see Fig. |l|, left plot). As previously noted (see 
Sect. Q), the emission of YZ CMi is often considerably po- 
larized. Throughout this observation the polarization was 
predominantly left circular (about 60 % circular polarisa- 
tion during the quiescent emission, and 90 % during the 
flares) . Two strong flares appear in the lightcurve obtained 
with the VLA. Their flux values reach up to 6.0 mJy. The 
detection of YZ CMi with both instruments was clear, as 
the rms noise in the VLA image was 0.063 mJy/beam, and 
in our high sensitivity VLBI maps made using only VLBA- 
VLA baselines the noise was 0.13 mJy/beam. Given these 
map noise values the peak intensity of the star was 47 a 
in the VLA image, and 16 a in the VLBA image. 

The mean flux density value for AD Leo was at a more 
typical level of 0.7 mJy (see Fig. [jj right plot). Simi- 
lar fluxes at 18 cm have previously been reported (e.g. 
[Jackson et al. 1989 ; Bcnz et al. 1995). A weak flare ap- 
pears in the VLA lightcurve (2.1 mJy). The star was de- 
tected by both instruments despite its low flux. The rms 
noise for the VLA map was 0.035 mJy/beam, while for 
the map made with only the VLA- VLBA baselines it was 
0.18 mJy/beam. Peak map values were 15.6 a for the 
VLA and 5 a for the VLBA maps respectively. Tracing 
the lightcurve of AD Leo from the interferometric VLA 
data was a delicate task because of two strong sources in 
the field of view. 



3.2. Size and shape of YZ CMi 

The phase referenced VLBA map of YZ CMi is shown in 
Fig. H. Evidence for spatial resolution was found in this 
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Fig. 1. Total flux density (stokes I) against time during the YZ CMi and AD Leo observations, extracted from the 
VLA data. The relatively large average flux of YZ CMi is noticeable as are several flares. The horizontal numbered bars 
in the plot show the periods (scans) during which the VBLI data were collected. Each scan is composed of subscans 
on the target star and the calibrator alternately. 



map. Further evidence for extended emission can be most 
clearly seen by examing fringe amplitude versus baseline 
length obtained using the VLA-VLBA baselines (see Fig. 
||, left plot). We obtained this plot by first coherently av- 
eraging the data on each baseline-scan (75 minutes) in or- 
der to increase the signal to noise ratio significantly above 
unity. We then binned the amplitude over baseline length 
finding the mean amplitude in each bin by incoherent av- 
eraging. The error bars on this average were determined 
from the internal scatter of the data. 

The fall off noticeable in Fig. || left plot clearly indi- 
cates a resolved source. What is more, the phase values 
on VLBI baselines to the phased VLA over the whole ob- 
servation show no significant variation from zero. There 
is therefore no evidence for anything other than a single 
centro-symmetric component. We searched also for other 
evidence for non-symmetrical structure looking at closure 
phases, but the SNR of these were too low. 

Given the close to quadratic fall off of the fringe am- 
plitude with M,w-distance shown in Fig. || left plot it is 
impossible to distinguish between gaussian, sphere, disk 
or ring like models (Pearson 1995). We therefore fitted 
one-component gaussian models to the YZ CMi data. The 
numerical values for the fits are summarized in Table |l|. 
The dimensions of sphere, disk and ring which would show 
similar fits are expected to be respectively 1.8, 1.6 and 1.1 
times the gaussian FWHP values (Pearson 1995). Two 
ways of fitting were followed: fitting in AIPS using the 
task UVFIT and our own model fitting to the data. The 
first fitted an elliptical gaussian and obtained for the whole 
data set a major axis of FWHP of 1.4 ±0.3 mas and a mi- 
nor axis of 0.5 ±0.25 mas (la errors). There is therefore 
no strong evidence for ellipticity, and our subsequent fit- 
ting of the data outside of AIPS assumed only a circular 



gaussian. With such a model it was possible to fit most of 
the data within 1 a. The best fitting FWHP size of the 
corona was found to be 0.98 ±0.2 mas, which corresponds 
to 1.7 ±0.3 stellar diameters. 

Since the second VLBI scan corresponds precisely to 
one of the two strong flares (see Fig. |l|, left plot), it was 
interesting to study it more closely. Fig. || (right plot) 
shows the amplitude versus u,v-distance. Of this scan we 
selected the subscan on the target with the highest flux 
density value and coherently averaged the data over it (3 
minutes), obtaining one point per VLA-VLBA baseline. 
The solid line corresponds to the same gaussian model 
fitted to the whole data set (see Fig. ||, left plot) except 
for the total flux density which was increased to 5.1 mJy 
We find that this scaled model fits the data within the 
errors and therefore there is no evidence for a change in 
source size during the flare. Independent gaussian fits to 
the data are also consistent with this conclusion (see Table 

0). 

We should add that the contribution of the proper mo- 
tion and of the changing parallax of the star during the 
ten hours of observation is -0.32 mas and -0.16 mas in a 
and 5, respectively. These values are small enough that 
they do not give a significant contribution on the spatial 
extent. 



3.3. AD Leo 

An image of AD Leo is shown in Fig. [| It appears 
slightly elongated. This image is probably affected by the 
star's high proper motion (Table [l]), since the extension 
is exactly along the expected direction: the star appears 
blurred because of its motion during the synthesis obser- 
vation. 
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Fig. 2. A contour plot of YZ CMi observed on April 18/19, 1997 (left) and AD Leo observed on December 12, 1996 
(right). The maps have rms noise cr=0.13 mJy/beam and a— 0.18 mJy/beam respectively. The first contour is at 3 a, 
and then at steps of one a. At the bottom left of each panel the clean beam is drawn, and the circle in the lower right 
corner represents the expected optical size of the star, indicated in Table ^. 




Fig. 3. Fringe amplitude versus u,u-distance of YZ CMi on the sensitive VLA-VLBA baselines. The solid line is the 
result of fitting a circular gaussian (see also Table Q). Note the change in vertical scale of the right plot compared to 
the left one. The right pannel shows the amplitude versus u, v-distance during the strong flare in scan 2. The solid 
line in this plot represents the same gaussian model used in the left plot except for the total flux density value which 
has been increased to 5.1 mJy 



As shown in the flux curve of AD Leo (Fig. right 
plot), the total flux varied by a factor of 3 during the ob- 
servations. Adding this fact to the high proper motion of 
the star, it is not surprising that attempts to fit a single 
gaussian to the u,u-data did not give consistent results. 
The weakness and variability of the star make the estima- 
tion of errors on the size difficult. However, from the image 



(see Fig. ||), we can note that the apparent FWHP per- 
pendicular to the motion corresponds to the beam FWHP 
in this direction. The intrinsic FWHP size of the emitting 
region is therefore likely to be less than half the beam 
FWHP or about 1 mas, which equals the estimated opti- 
cal diameter of the star (see Table |l|) . This might indicate 
a very compact corona or an emitting spot on the surface 
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YZ CMi 1997 


AD Leo 1996 


Spectral type 


M 4.5Ve 


M 4.5Ve 


Distance 


5.93 ±0.09 pc 


4.90 ±0.07 pc 


Optical diameter 


0.58 mas 


0.95 mas 


(UVFIT all obs) 


1.4x0.5 (±0.3x0.25) mas 






PA 12.8° (±5.1) 




(gaussfit all obs) 


0.98 (±0.2) mas 




(UVFIT scan 2) 


0.890 (±0.22) mas 




(gaussfit scan 2) 


1.02 (±0.2) mas 




pm and plx 


Aa: -0.32 mas 


Aa: -0.90 mas 


during observations 


AS: -0.16 mas 


A<5: ±0.29 mas 


VLBA Position 


a: 7" 44™ 40!22399 ±0.00007 


a: 10 h 19 m 36=40112 ±0.00007 


(J2000) 


5: 3° 33' Iff! 1140 ±0.0005 


5: 19°52'12'.'1008 ±0.0005 


Mean flux density 






during observation 


2.97 mJy 


0.54 mJy 


Brightness temperature 


7.29-10 7 K 


> 4.93 • 10 7 K 



Table 1. Summary of the results from the observations and fits for YZ CMi and AD Leo. The sizes of the radio 
emission are labeled with the name of the procedure used to obtain them (UVFIT in AIPS, gaussfit outside AIPS). 
The distances are taken from the Hipparcos and Gliese catalogues, respectively, and the optical diameters are from 
Pettersen (1980). For the AIPS fit of scan 2 we constrained the program to fit a circular gaussian because of the small 
number of visibilities. We also did a fit for AD Leo but it did not converge (see text). 



Star 



Name of calibrator 



*(J2000) 



<5(J2000) 



YZ CMi 



0763±017 
0743-006 



07 h 39" 
07 h 45" 



18.0339 s 
54.0823 s 



01° 37 04.619 
-00° 44 17.538" 



AD Leo 



1022±194 
1013±208 



44.8096 s 
44.3251 s 



19° 12 20.416 ' 
20° 37 47.336" 



Table 2. Adopted positions for the correlation of the primary and secondary calibrators of each target star. The 
primary calibrators were used as reference sources in the phase referencing, the secondary calibrators as check sources 
(sec Sect. 2). All positions come from Johnston et al. (1995). 



of the star. A very conservative upper limit on the size 
of the corona in the former case would be to assume the 
emitting region to be an optically thin sphere instead of a 
gaussian (Pearson 1995) in which case the diameter is less 
than 1.8 times the photosphere diameter, and it therefore 
has an extent above the photosphere of less than 0.8 i?*. 



3.4- Astrometry 

The relatively large signal-to-noise ratio for YZ CMi and 
the phase referencing to a calibrator with good (<0.5 mas) 



positional accuracy in the radio frame ( Johnston et al 



1995) allowed us to determine a precise position for this 
star. This position was compared with the position given 
by the Hipparcos catalogue (ESA 1997). Correcting for 
proper motion and parallax, we found a discrepancy of 
20.9 mas in a and 30.4 mas in 5, thus a total deviation 
of 36.9 mas. The proper motion of the star is given in 
the Hipparcos catalogue as -344.9 ±2.6 mas/yr in a, and 
-450.8 ±1.75 mas/yr in 5. Considering the time interval 



between the two measurements of 6 years, the difference 
is 2 a and thus within the accuracy of the proper motion 
error bars. Combining the VLBA and Hipparcos positions 
(courtesy of F. Arenou), an improved proper motion of 
-348.6 ±0.6 mas/yr in a, and -446.6 ±0.3 mas/yr in S can 
be derived. 

The position of AD Leo obtained with the VLBA was 
compared with those available in the Gliese and Tycho 
catalogues. The latter showed a deviation with the VLBA 
position of 176.3 mas and 100.0 mas in a and <5, respec- 
tively. They are within one standard deviation of the Ty- 
cho catalogue accuracy. 

4. Discussion and conclusions 

VLBA observations have spatially resolved YZ CMi and 
the data could be fitted with a circular gaussian of a 
FWHP of 0.98 ±0.2 mas. The radio corona extent is 
1.77 x 10 10 ±8.8 x 10 9 cm above the photosphere (the pho- 



spheric radius is assumed to be 2.6 x 10 cm, Pettersen 



1980). For AD Leo, which is closer and has a larger photo- 
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star radio luminosity extent of corona above photosphere 

[erg/s Hz] cm R+ 



reference 



Sun 

AD Leo 
UV Cet B 
YZ CMi 



5.7 x 10 iu 
3.04 x 10 13 

8.8 x 10 13 
1.4 x 10 14 



7.4 x 10 s 
< 2.8 x 10 1 
2.1 x 10 10 
1.77 x 10 10 



0.011 R+ 
< 0.8R+ 
2.1 R* 
0.7 R* 



Aschwanden et al 



1995 



this work 



Benz et al 



1998 



this work 



Table 3. Luminosity and coronal size of Sun and dMe stars at 3.6 cm wavelength. 



sphere, but was observed at a much weaker flux level, the 
corona was not resolved and we set a robust upper limit 
of 2.8 x 10 10 cm above the photosphere (see Sect. 3.3). 

4--1- The brightness temperature 

For YZ CMi (see Table Q) we obtain a mean bright- 
ness temperature Tf, = 7.3 x 10 7 K, while for AD Leo 
we set an upper limit of T b = 4.93 x 10 7 K. These 
mean Tf, values are smaller than previously reported for 
YZ CMi by [Benz fc Alej (1991) and AD Leo at 18 cm 
( Jackson et al. 1989J ). The lower values at 3.6 cm are still 
consistent with a non-thermal spectrum from gyrosyn- 
chrotron but do not formally exclude thermal processes. 
However the significant circular polarisation found dur- 
ing the observations strongly argues for a gyrosychrotron 
emission mechanism. 

The derived extent of the coronae above the photo- 
sphere of the dMe stars is compared to the Sun in Table 
|3|. The solar value refers to stereoscopic measurements of 
the thermal gyroresonance emission of active regions. The 



average value at 10-14 GHz reported by |Aschwanden et al. 
(1995) has been used. 

The experience from the solar radio emission makes it 
clear that the observed radio size is only a lower limit of 
the size of the stellar corona. Nevertheless, these results 
indicate that the observed active dMe stars have much 
larger active coronae than the Sun. It might be that such 
dMe stars posess systems of closed loops reaching heights 
in excess of a stellar diameter. One way to realize such 
extended coronae is by large distances between footpoints 
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as possibly seen in the case of UV Cet B (Benz et al. 1998). 
This indicates either that active regions are very large, 
or that active loops preferentially connect different active 
regions. 
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